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Abstract 

We consider the Navier-Stokes system describing the motion of a compressible barotropic fluid 
driven by stochastic external forces. Our approach is semi-deterministic, based on solving the system 
for each fixed representative of the random variable and applying an abstract result on measurability 
of multi- valued maps. The class of admissible driving forces includes the (temporal) white noise and 
the random kicks, considered recently in the context of incompressible fluid models. 
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1 Problem formulation 

We consider the Navier-Stokes system governing the time evolution of the density g and the velocity u 
of a compressible viscous fluid, driven by a stochastic external force that can be formally written in the 
form: 

dg + div x (gu) dt = 0, (1.1) 

d(pu) + (div x (gu ® u) + V ' x p{o) — div^S^ dt = gdvs. (1-2) 

The quantities g = g(t 1 x 1 uj) 1 u = u(i, x, u>) are functions of the time t € (0, T), the spatial position 
x € Q, and to E O = {O, B, /.i}, where O is a topological probability space, with the family of Borel sets 
B, and a regular probability measure fi. System (jl.ip . (|1.2p is supplemented with the standard no-slip 
boundary condition 

u| an = 0, (1.3) 
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and the initial conditions 



0(0,-) = Qo, gu(O,0 = (£>u) Q . 



(1.4) 



The symbol § denotes the viscous stress determined by Newton's rheological law 

§ = ^(V^u + V^,u - ^div^u) + TjdivaUl, f > 0, rj > 0, (1.5) 

p = p{g) is the pressure, and the perturbation w is a random variable represented for a. a. w by a bounded 
function, sufficiently regular with respect to the spatial variable x £ fl. 

There is already a substantial amount of work on the incompressible stochastic Navier-Stokes sys- 
tem, see the surveys by Bensoussan [2j, Mattingly [21], and the references cited therein. In 3D general 
existence results for equations with Gaussian noise have been obtained for the so-called weak solutions 
in probabilistic sense (or martingale solutions), cf. [13], [13] or [IS]. Uniqueness of solutions in 3D is 
open just as in the deterministic case, but working carefully with the so-called Markov selections ([7|, 
[15) ) it was possible to prove existence, uniqueness, ergodicity and strong mixing of the invariant measure 
(stationary solution) assuming certain nondegeneracy of the noise (for a general result on Markov selec- 
tions including also the 3D stochastic Navier-Stokes equation cf. [E]). Analogous results (existence of 
martingale solutions and Markov selection) for 3D Navier-Stokes equations with jumps (driven by Levy 
noise) have been recently proved by Dong and Zhai [9]. 

Much less seems to be known for compressible fluid flows. To the best of our knowledge, the only 
available results concerning the compressible stochastic Navier-Stokes system concern the 1-D case, see 
Tornatore , Fujita Yashima [24], [25], and the rather special peridic 2-D case examined by Tornatore [23] 
by means of the existence theory developed by Vaigant and Kazhikhov |26] . 

Similarly to the seminal paper by Bensoussan and Temam [3], our approach is semi-deterministic 
based on the concept of weak solutions to the Navier-Stokes system, in the framework of the existence 
theory developed by Lions [2U] , with the extension of 1TJ . The class of the weak solutions is sufficiently 
robust to incorporate the driving forces f of very low regularity with respect to the time, in particular, 
we may consider 

f = f (t, x) = dtvr(t, x), as soon as w is Lipschitz with respect to the spatial variable x € f2. 

Thus the stochastic problem may be solved "pathwise" for any individual choice of the random element 
u>. It is also shown that under some mild restrictions on the paths of driving process ( the time integral 
of the noise) the solution is a "stochastic process", i.e. it is a (measurable) random variable with values 
in the paths space. As an example we consider the Levy noise , which is important because it represents 
a general stationary incorrelated noise driving the paths which are stochastically continuous (cf. Section 
5 for details). 

1.1 Weak formulation 

Taking advantage of the specific form of (|1.1[) . (I1.2[) . we may formally write 

d(f>u) — gdw = d[g(u — w)] + dgw = d[g(u — w)] — div K (giu)w dt. 

Accordingly, the product div x (gu) may be interpreted in the sense of distributions provided w ranges in 
the Sobolev space Wj' 9 (fi; i? 3 ) for a sufficiently large exponent q. 
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We say that a pair of random variables g, u, 

g = g(t,x,uj), Q (;0j)€C([0,T];L\n)), (1.6) 

u = u(t,x,u), u(-,w) € i 2 (0,T; W^ 2 (ri;R 3 )) for a.a. u G O, (1.7) 
is a weafc solution to the problem (|l.ll - [L^]) , if the following relations hold: 

Equation of continuity: 

1. 

Q > 0, p£ L°°(0,T : L 7 (ft)) for a certain 7 > 1 for a.a. u) G O, (1.8) 
and the family of integral identities 

J (^q(t,x,lj) + b(g(T,x,u})) S Jip(x) dx - J (^Qo{x, uj) + b(g (x, oj)fj ip(x) dx (1.9) 

— J J (^g(t,x,u)u(t,x,u))-V x p(%)+(^b(g(t,x,uj))—b\g(t,x,(d))g(t 

is satisfied for any test function ip G C 1 (ri), any b G C^°[0, 00), and a.a. ui G 0; 

Momentum equation: 

2. 

G L°°(0, T : L 2 ^ +1 )(f>; i? 3 )), g(u - w) G C7 woak ([0, T]; £ 2 ^ +1 )(fi; i? 3 )) for a.a w G O, 

(1.10) 

G L°°(0, T; L 1 ^) n L 9 ((0, T) x n) for some for a.a. uj G O, (1.11) 
and the integral identity 

g(r, x, uj) I u(t, x, uj) — w(r, a;, uj) j • dx — / (gu)o(a;, w) ■ <p(x) dx (1-12) 

g(t, x, cu)u(t, x, uj) ® u(i, a;, w) : V x tp(x) + p(g(t, x, uj))div x p(x)) dx dt 

dx dt 



JQ 



(§(Va;U(i, x,uj)) : V x ip(x) + g(t, x, w)u(t, x, uj) ■ V x (vr(t, x, uj)ip(x) 



'0 JQ 

?3 



holds for all cp G C c °°(fi; R 6 ), and for a.a. uj G C 

Energy inequality 

3. 



^g(r,a;,a;)|u(T,x,£j) - w(r,x,w)| 2 + P(q(t,x,oj)) ) d.r (1.:13) 
§(V x ii(i, x, a;)) : V s u(t, x, w) dx 



n 



/ ( ^q(s,x,uj)\u(s,x,uj) - w(s,x,uj)\ 2 + P(g(s,x,u))) j d.r 
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+ Is I i^^ xuS> ' VxW ~ gu ® u : ^ xW ~ p(e)^ v xw + ^gu ■ V x \w\ 2 ^ (t, x,u) dx dt 
for a. a. r > s > including s = 0, and a. a. ueO, with 

P(Q) = Q f dz. (1.14) 

Remark 11.11 1 The family of integral identities fl.9\) represents a weak formulation of the renor- 
malized equation of continuity introduced by DiPerna and Lions JEI- The more standard, but in fact 
equivalent, "distributional" formulation of il.9\) . ll.lSfy . U.13\) will be given in Section^ below. 

We can check easily that for the energy inequality (j!.13[) to yield uniform bounds on the family of 
solutions, it is necessary that 

ess sup ||w(i, -,uj)\\ w i,^ (Q . Ra) < c(uj) for a.a. uj e O, (1-15) 

te(o,T) 

in particular, certain regularity with respect to the spatial variable is needed. 



1.2 Main results 

Our main goal is to develop an existence theory for the Navier-Stokes system (|1.8I - 11.13[) . where the 
solutions q(-,lj), u(-,w) are random variables ranging in the function spaces specified through (jl.6|) . 
(11.71) . Analogously to Bensoussan and Temam [3], our approach is based on: 

• a general stability result for solutions of (|1.8I - H~T3]) . where the forcing term satisfies only p. 151) ; 

• application of the abstract measurability theorem for multivalued maps, similar to Bensoussan and 
Temam [3J Theorem 3.1]. 

By stability we mean that a sequence of weak solutions to (|1.8I - [L13]) . with pre-compact data go, 
(f?u)o, w, admits a subsequence that converges (weakly) to another solution of the same problem. Besides 
the nowadays standard ingredients of the existence theory (see [TT], Lions [2D]), the proof is based on 
compactness of certain commutators in the spirit of Coifman and Meyer [B] . 

The paper is organized as follows. In Section [5] we recall the standard definition of weak solutions to 
the Navier-Stokes system and prove a general stability results assuming low regularity of the driving force 
represented by w. This result is exploited in Section [3j where we show existence of the weak solutions 
for system (|1.8I - H~T3]) . with a fixed (irregular) force w. Section |4] is devoted to the analysis of the 
associated stochastic system. We introduce the suitable function-spaces framework and use the abstract 
measurability theorem to prove that the weak solution generate a random variable. Several concrete 
examples of the driving force w are discussed in Section [S] 

2 Weak sequential stability for irregular forcing 

In order to fix ideas, we suppose that the pressure p belongs to the class C[0, oo) D C 2 (0, oo) and satisfies 
p(0) = 0, p'(g) > for g > 0, lim ^4 = P~ > 0, 7 > 3/2. (2.1) 

e^oo gl 1 

In this section, we drop the parameter uj and fix a function w. Accordingly, we may rewrite the weak 
formulation of the problems as follows: 
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Equation of continuity (renormalized): 
[ [ ((0 + 6(e)) %>+ (e + &(*?)) u-V^) dx dt (2.2) 
(b'(e)e - di Va; u) 93 dx dt - / (g + 6( eo )) ^(0, •) 



d.r 



/o Jo 

for any test function tp e C£°([0, T) x H), and any 6 G [0, 00); 

Momentum equation: 

2. 



y y (g(u — w) • d t f + gu (g) u : V x </3 + p(g)div x ( / 9^ da; di (2.3) 

/ / (S(V K u) : V x tp(x) + gu ■ V x (w ■ tp) J dx dt - / (gu) ■ ip(0, ■) dx 
Jo Jo v 'Jo 



for all ip e C7 c °°([0,r) x Sl;R 3 ); 

Energy inequality: 

3. 



f f (lg\u-w\ 2 + P(g)) dxd t ipdt+ [ [ S(V x u) : V x u da; ^ dt 
Jo Jo / Jo Jo 



(2.4) 



+ J j (§(V x u) : V^w- guigiu : V x w - p^div^w + ^gu ■ \7 x \w\ 2 ^j dx ip dt 

for any tj) £ C* c °°[0,r), V > 0. 

The heart of the paper is the following result that may be of independent interest. 

Proposition 2.1 Let SI C R 3 be a bounded Lipschitz domain. Suppose that the pressure p satisfies \2.1]) . 
Let {w„}" =1 be a sequence of functions, 

ess sup \\^n(t,-)\\ w i^ (n . R3) <C w ,w n ^w l nL 1 {0 7 T;W 1 ' 1 {Sl;R 3 )). (2.5) 

te(o,T) 

Let {gnjUnj^L-L be a sequence of weak solutions of the Navier-Stokes system in (0,T) x Sl 7 driven by w„, 
and emanating from the initial data £>o,n> iff u )o,n such that 

go.n > 0, / g , n dx = M > 0, ||go,„|U-r ( n) < %o, f l(gu)o ' ra| dx < E , (2.6) 
Jo Jo Bo,e 

g , n -> g in IS (SI), (gu) ,„ -> (Ho wea% m L 1 ^), / l(gu)o '" 12 da; -> / K^li^ d.T. (2.7) 



O 00," JO 00 
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Then, at least for suitable subsequences, 

Qn -> gin C weak ([0,T];iT(O)) and m L\(0,T) x fi), (2.8) 

u„ -> u iweoWy in i 2 (0, T; Wq' 2 ^; i? 3 )), (2.9) 

g n (u n -w„)^ e (u-w) mCweakC^T];^/^ 1 )^;^ 3 )), (2.10) 



where g, u is a weak solution of the same problem, with the driving force w, and the initial data go, 
(gu) . 

Remark [2jl In hypothesis V2.6\) , we tacitly assume that (gu)o, n {x) — on the set where go,n{x) = 0. 
The remaining part of this section is devoted to the proof of Proposition 12.11 

2.1 Uniform bounds 

Since the total mass M is a constant of motion, we get 

||0n(V)IU*(n) =M for any t G [Q,T]. (2.11) 
Consequently, in accordance with hypothesis (|2.5I) . 

ess sup / £> n |w„| 2 da; < c(M, (?„,); 
te(o,T) 

whence the standard Gronwall argument can be used in the energy inequality (|2.4[) to obtain 

ess sup \\yfg^u n (t,-)\\ L 2 (n . R3) < c{M,C w ,E Q ,T), (2.12) 
te(o,T) 

ess sup \\g n (t,-)\\ LH n)<c(M,C w ,E ,T), (2.13) 
te(o,T) 

and, by virtue of Korn's inequality, 

/ IMV)||^,= (n . fl3) < c(M,C w ,E ,T). (2.14) 

«/ 

The estimates (|2.11I - [2~T4| are exactly the same as the a priori bounds available for the compressible 
Navier-Stokes system, cf. |llj . In particular, we may infer, at least for suitable subsequences, that 

g n -+ e m C wcak ([0, T]; IS (SI)), (2.15) 

u„ -> u weakly in L 2 (0,T; Wq' 2 (SI] R 3 )), (2.16) 

and 

g n u„ -> wcakly-(*) in L°°(0, T; i 27 ^ 7+1 ^(r2; i? 3 ). (2.17) 
Moreover, a short inspection of the momentum equation (|2.3j) yields 

g„(u„ - w„) -> g(u - w) in C wca k([0, T]; I^^ +1 \S1; R 3 )). (2.18) 
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Finally, combining (I2.17[) , (|2.18l) with hypothesis (|2.5[) , we obtain that 

{g n u n ){t, •) -> (gu)(t, •) weakly in L 2rf/ ^ +l \Q; R 3 ) for a.a. t € (0,T); (2.19) 
in particular, as — 2j > |, 

g n u n -> gu in L q (0, T; W~ 1,2 (Tl', R 3 )) for any 1 < q < oo. (2.20) 
Relations (|2TTd| . (j2"T2T))) . together with the bounds established in (j2~T2"l) . (|2~T5|) . give rise to 

£>„u„ ®u„ pu (g) u weakly in L 2 (0, T; L 67/(47+3) (ft; i? 3x3 )). (2.21) 

2.2 Pressure estimates 

In view of the uniform bounds established in (|2.11l - l2.14p . the pressure p can be estimated exactly as in 
[12] . specifically, by means of the quantities 



Ken) - T^T J b(0n) dx 



i>eC™(o,T), b(e)*e u 



used as test functions in the momentum equation (|2.3[) . where B ~ div -1 is the so-called Bogovskii 
operator, see Bogovskii 5 , Galdi |16j . Indeed the extra terms that appear in (|2.3p in comparison with 
the situation treated in [12] are regular as a consequence of hypothesis (|2.5|) . Consequently, exactly as in 
[12], we deduce that 

T 

p(g n )gi dx < c(M,C w ,E ,T) for a certain (3 = f3(M,C w ,E ) > 0. (2.22) 



o Jn 



2.3 Pointwise convergence of the densities 

In order to complete the proof of Proposition ^. 11 we establish the following crucial result: 

Q n -> g a.a. in (0, T) x Cl. (2.23) 

Obviously, the a.a. pointwise convergence of the densities is a key ingredient of the existence theory for 
the Navier-Stokes system, see [H], Lions [20]. Here, the situation is even more delicate because of the 
presence of the new terms containing w. 
Following the arguments of [TT], we use 

<p = ^MaOV.A-^ft,)], il> £ C c °°(0,T), r,,£ e C c °°(«) 

as test functions in the momentum equation (|2.3[) . and, similarly, we take 

^ = V^)V. A" 1 [£&(>)] 

in the limit equation 

J J (q(u — w) • d t cp + gu u : V x <p + p(g)div x (p^ dx dt (2.24) 
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(w • ip) dx dt — / (gu)o ■ ip{0, ■) dx, 



= [ [ §(V K u) : V x ip dx dt+ [ [ gu ■ V a 
Jo Jn Jo Jo. 

where we have used h(g) to denote a weak limit of the compositions h(g n ). 
Following step by step the arguments of [IT] we deduce that 

lim / / ^Up{g n )b{g n ) - S e : V ^A^V x [^b(g n )]) dx dt 

i-n(W)W)-^ ■ Vx^-^x^Ko)}) dt 

• rjg n u n — £&(g„)V x A V x [qg n u n ] ) dx dt 
l ^°°Jo Jr. 3 v y 

-jf y ^u- (VsA^VsK&te)] ■i/ju-^V I A- 1 V ai [f)gu]) dx dt 
^^nUn • V x (r/w„ ■ V^A -1 [£&(£„)]) + g„w„ • <9 t (^V^A -1 [£i>(e n )] 



(2.25) 



T 

./ii 3 



lim 



Ji? 3 
T 



dx dt 



o ^i? 3 



da; dt. 



2.3.1 Weak stability of the effective viscous flux 

We show that the right-hand side of (|2.25l) vanishes, in particular, we recover the so-called weak continuity 
of the effective viscous flux, first established by Lions [20 : 

p{e)b(g) - + Xj b{g)div x u =p(ojKQj-(^v + x\ b(g)div x u (2.26) 

for any bounded function b. In order to see (|2.26|) . observe first that we can use the same arguments as 
in [TT] to show that 

lim / / Vu„ • ( V x A^ 1 \7 x [^b(g n )} ■ T]g n u n - £b(g n )V x&^V x [r]g n \i n } ) dx dt (2.27) 

n^ocJ Q J R3 \ J 



T 

JR S 



i^u ■ (v x A- x V x [e&(e)] • VP" ~ eK^VzA^V^wu]) dx dt. 

Indeed the only difference between |llj and the present situation is that the convergence 

g n u n -> gu in C wcak ([0, T]; L^ +1 \Sl)), 

valid in [11], is replaced by a weaker statement (|2.19|) . Fortunately, relation (|2.19p is still sufficient to 
carry over the proof of (|2 .27[) in the same way as in QT] . 
Thus our task is to show that 



lim 

n— ¥oo 



'0 J R 3 



?Pg n u n ■ V x (r]w n ■ V X A ^(fc,)]) + g„w„ • d t (^VxA -1 dx dt (2.28) 
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da; dt, 



da; dt 



(2.29) 



da; dt. 



~i>QU ■ vjrfw • VsA- 1 ^)]) + gw ■ ftf^V^A" 1 ^)] 

/O J_R3 L V ' V 

which, after a straightforward manipulation, reduces to 

lim ipr) p„w„ • (u„ • V^A^divz [£&(£„)] - \7 x A~ 1 div x {£b(Q n )u n ] 

n ^°°Jo Jr3 L V 

= / jL» ^ ' ( U ' V * A ~ ldiv *K^>] - VxA^divxK^u]) 
Now, with (|2~T4|) at hand, we observe that the commutator 

(u„ • V x A _1 div x [£&(£ n )] - V a A~Miv x [£&(£„)u n ]J 

is actually bounded in L 2 (0, T; W 1,q (R 3 ; R 3 )) for any 1 < q < 2. This is easy to see computing div and 
curl of this quantity and using continuity of the Riesz operator in L q — spaces, see also the first part of 
the proof of Theorem 10.28 in [ID]. 

On the other hand, as a direct consequence of (I2.5[) . (|2.15p : 

g n w n -> in L 2 (0,T; W^'(0;i? 3 )); 
whence ((2~2"D|) follows. We have proved (|2~2l)|) . 



2.4 Conclusion 



Having established (|1.8j) we can finish the proof of (|2.23j) by the same arguments as in [IT]. Moreover, 
the pointwise convergence of the densities may be strengthened to 



g in C([0, T];L q (£l)) for any 1 < q < 7. 



(2.30) 



In particular, p(g) = p(g). We have proved Proposition 12.1 



3 Existence of solutions for problems with irregular forces 

Our aim is to use Proposition 12 . 1 1 to show existence of weak solutions to the Navier-Stokes system driven 
by w, 

ess sup ||w(t,')|| w -i,oo( n . ii 3) < C w . (3.1) 
te(o,T) 

To this end, suppose first that w is regular, more specifically, 

weC'dO.r];^!!;^), w(0, •) = 0. (3.2) 

As shown in [TT] Theorem 1.1], the Navier-Stokes system admits a weak solution in the sense specified 
at the beginning of Section (3] where the momentum equation (|2.3p is replaced by the integral identity 



J J (^gu ■ d t f + gu <g> u : V x i/; + p(g)div x ipj dx dt 



(3.3) 



!) 



T 

JQ 



(s(y x u) : V x p(x) - gd t w ■ tpj dx dt - / (gu) ■ p(0, •) da; 
for all <p € C%°([0,T) x 17; i? 3 ), while the corresponding energy inequality reads 



/ (^g\u\ 2 + P(g)) da; dt + [ [ §(V*u) : V K u da; V dt (3.4) 

JQ \* / JO Jf2 

<^(0) / f^ l(gu) ° |2 +P( g0 )) dx+ [ T [ eftwudarVdt 

for any r/> g C£°[0,T), ip > 0. The result, as stated in [TT], requires certain smoothness of the boundary 
<9£1 that can be relaxed, see Kukucka [T5]. In particular, the existence can be shown provided ft is a 
bounded Lipschitz domain. 
By virtue of the identity 

/ / gd t w ■ p dx dt = / / gd t (w ■ tp) dx dt — / / gw ■ d t <p dx dt 
Jo JQ Jo Jn Jo Jq 

= — J gu- V x (w ■ p) dx dt — / / gw -dtp dx dt 
Jo Jn Jo Jn 

satisfied for any tp € C^°([0, T) x Ct; i? 3 ), relation J33]) reduces to 
Moreover, taking <p = w as a test function in (|3.3p yields 



(3.5) 



y y ^pu • <9fW + gu ® u : V K w + p(g)div x w S J dx dt 

= 11 (§(V a u) : V K w - gi9 t |w| 2 ) da; dt + / (gu)(r, •) • w(r, •) da; 
Jo Jn K A i J n 

J f Q^lwl 2 - Q u • w ) ( T , ') da; + J J Qgu • V x |w| 2 + E>(\7 x u) : V a wJ da; dt for any r > 0, 



which, together with (|3.4|) . gives rise to the energy inequality (|2.4p . 

Having obtained solutions for any regular w we may approximate a general function w by a sequence 
of smooth functions w n satisfying (|2.5[) and use Proposition 12. II to conclude: 

Proposition 3.1 Let O be a bounded Lipschitz domain. Suppose that the pressure p satisfies i2.1\) , 

weCaktO^;^' 00 ^^ 3 )), (3.6) 

and 2 

g Q e L 7 (^), go > 0, / g dx = M > 0, {gu) E L\n;R 3 ), [ ^ 0U ^ dx < oo. (3.7) 

JS1 Jf2 f?0 

TTien t/ie Navier-Stokes system admits a weak solution in (0, T) x in the sense specified in \2.2\ - 
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Remark [3jl The symbol L"^ cak denotes the space of weakly measurable functions. 
Remark [3j 2 Writing 

(£»U) = y/Q0- 



we can see that \3. 7| ) implies that 



(fm) eL 2 ^ +1 >(ft;i? 3 ). 



4 Problems driven by stochastic forces 

We start by identifying the function space VV for w. We suppose that 
1. W is a separable complete metric space; 
2. 

WcC ak (0,T;W ^(l];i? 3 )); (4.1) 

3. if 

w n ->win W, 

then 

ess sup ||w n (t,-)|| w i.« m . fl 3, < c, and u>„ -> u; weakly - (*) in L°°(0,T; W 1,0o (f2)). (4.2) 
te(o,T) ^ 

Next, we introduce the spaces of the data. In accordance with hypotheses of Proposition 13. 11 we take 

td = {[eo,(eu)o]| eo e £ 7 («), (eu) e L 2 ^ +1) (^;i? 3 ), (4.3) 

go > 0, / go dz = M > 0, / ^ l(gu) ° |2 + P(q ) dx < E Q } . 
Jn Jn * Qo ) 

We note that XV is a closed convex subset of the separable Banach space L 7 (f2) x L 2l ^ J+1 ^ (Q; i? 3 ), in 
particular, it is a Suslin space. 
We suppose that the mapping 

lj E O ^ [g (-,oj),(gu) {-,u)} elV 

is a random variable on O = {0,B,fj,}, with a regular probability measure [i. Thus for any e > 0, there 
exists a compact set K, e C O such that 

(i(0\lC e )<e, 

and the mapping 

cj h- > [go( - ,w), (gu)o(-, cj)] restricted to K, e is continuous. 
Similarly, we assume that 

weOn> w(-,lj) e w 
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is a random variable. 

Finally, we consider a multi- valued mapping 

S: {[f>o,(f?u)o],w} elVx [g, u ] G C([0, T]; i x (Q)) x Z 2 (0, T; W^O; -R 3 )) 

that assigns to the initial data £o> (£ ,u )o an d to the forcing w the associated family of weak solutions 
to the Navier-Stokes system, the existence of which is guaranteed by Proposition 13.11 Similarly to the 
situation examined by Bensoussan and Temam [3], the weak solution need not (is not known to be) be 
unique. As a consequence of Proposition the values of S are non-empty closed subsets of the Banach 
space C([0,T];i 1 (n)) x L 2 (0, T; W 1,2 (fi; i? 3 )), and S possesses a closed graph. 

Applying the abstract result of Bensoussan and Temam J3J Theorem 3.1, Lemma 3.1] we conclude that 
for each oj G O, there exists a weak solution of the Navier-Stokes system in the sense specified through 
(|1.8I - II . 12|) such that the mapping 

«-> [ e (.,w),u(-,«)] G C([0,T};L\n)) x L 2 (0, T; W 1,2 (f2; R 3 )) 

is a random variable. 

We have shown the main result of the present paper: 

Theorem 4.1 Let fl C i? 3 6e a bounded Lipschitz domain. Suppose that the pressure p = p(g) satisfies 
hypothesis R2.1\) , and that W is a complete separable metric space satisfying \4-<fy - Let 

ljeO^ {[eo(-,w),(gu)o(-,w)],w} elVxW 

be a random variable. 

Then there exists a random variable g(-,uj), u(-,w), 

q(-,w) G C([0,T];L\n)), u(-,w) G L 2 (0, T; W 1,2 (f2)) /or a.a. .eO, 

satisfying the Navier-Stokes system il.8\ - [LT3\) for a.a. u) G C 

5 Applications 

We conclude the paper by presenting a more specific type of random perturbation that satisfies the 
general assumptions stated in Section 21 which is the so-called Levy noise. For the convenience of the 
reader, the definition and basic properties of this type of random perturbations as well as some basic 
examples fitting into this framework are recalled. For more details see e.g. the basic monographs [22] or 

Let (0, J 7 , P) be a probability space and E — (E, || • a separable Banach space. 

Definition 5.1 An E-valued Levy process is a stochastic process L — (L(t),t >0),L: [0,oo) xO->£ 
satisfying: 

(i) L has independent increments, meaning for each < t\ < ti < . . . t n , n G N, the E-valued random 
variables L(t\) — L(to), L(t2) — L{t\), . . ., L(t n ) — L(t„_i) are stochastically independent. 

(ii) L has stationary increments, meaning the probability laws of L(t) — L{s), t>s, depend only on the 
difference t — s. 
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(iii) L(0) = and L is stochastically continuous, meaning for each r > and s € [0, oo), 

\imP[\\L(t) - L(s)\\ E > r] = 0. (5.1) 

t—>s 

Note that if the other conditions of Definition 15.11 are satisfied, (|5.1[) is equivalent to 

lim P[||£(f)|U > r] = (5.2) 

*->o+ 

for each r > 0. 

Recall that a process is called cadlag if its random paths are right-continuous and there exist left 
limits at each point. 

Theorem 5.1 Every E-valued Levy process L has a cadlag modification, that is, there exists a cadlag 
process L which is a Levy process itself (i.e. all random paths of L are [0,oo) — s- E cadlag) such that 

P[L(t) = Lit)] =1, te [0,oo). 

The proof is easily deduced from the classical Kinney's result on cadlag modification for Markov 
processes, cf. [15] . 

As is usual in probability theory we may identify a stochastic process with its modification, so with 
no loss of generality we assume in the sequel that L is cadlag. 

Example 5.1 E- valued Brownian motion with drift. Let L(t) =at + W[t), where a € E is fixed 
and W — (W(t),t € [0, oo)) is an E-valued centered Gaussian process with the incremental covariance 
operator R 6 C(E* , E), i.e. 

E(W(s),x) E , E *(W(t),y) E , E * = (Rx,y) E ,E>{sAt) 

for each x, y E E* , s,t € [0, oo). For example, if E is a Hilbert space, an operator R 6 £{E) is an 
incremental covariance of a Wiener process only if R is nonnegative, symmetric and trace class. The 
process L is then a Levy process (with continuous paths). It corresponds to the case when the equation is 
perturbed by noise white in time. 

Example 5.2 Let v be a finite measure on E \ {0}. An E-valued compound Poisson process with the 
Levy measure v is a cadlag Levy process satisfying 

oo , u 

P[i(()eA]= e -^» t ^-/(4), t>0, T£B{E)\{0}. 

k=0 

An equivalent definition that yields also the construction of the compound Poisson process is this: Let 
(l^) rie N be a sequence of independent E-valued random variables with identical distribution \v, b = 
v[E \ {0}), and assume that (Tl(t),t > 0) is a standard R + -valued Poisson process. Then 

n(t) 



Lit) =}Z k (5.3) 



k=0 



is a compound Poisson process with the Levy measure v. Clearly, L is a pure jump process with jump 
size Zk where II(t) is the number of jumps on the interval (0,t). If ~EL(t) — t J B >, yXv(dx) < oo, we 

may define the process L(t) :— L(t) — EL(i), which is called the compensated compound Poisson process. 
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The special Levy processes introduced in Examples 15.11 and 15.21 may be used to formulate a useful 
characterization of a general Levy process. 

Theorem 5.2 (Lcvy-Khinchine decomposition). Every E-valued Levy process has the representation 

oo 

L(t) =at+ W(t) + L n(t) + L (t), t > 0, (5.4) 

Tl=l 

where a G E, (W(t),t > 0) is an E-valued Brownian motion, L n for each n > 1 is a compensated com- 
pound Poisson process, Lq is a compound Poisson process, and W,Lq,L„ are stochastically independent. 
The series in (|5 ,4[) converges P -almost surely uniformly on each bounded interval o/R + . 

The processes Lq, L n may be defined so that their respective Levy measures are concentrated on the 
set {y; \\y\\ E > r }, {y; r n+1 < \\y\\ E < r n }, for any arbitrarily chosen sequence r„ \ 0. 

The natural state space for the paths on a time interval [0,T] of L is the linear space of E- valued 
cadlag functions T> — T>([0,T]; E). It may be equipped with the uniform norm but is not separable in 
this norm as needed in Section 01 A usual way to overcome this difficulty is to consider the so-called 
Skorokhod metric. Let us briefly recall its definition, for details we may refer to [H Chapter 3]. Let 
A denote the set of all strictly increasing, continuous mappings of (0,T] onto itself thus if A e A then 
A(0) = 0, A(T) = 1. Set 

d(x,y) = inf{e > 0| there exists A € A; X(t) —t\<e, \\x(t) - y(\(t))\\ E < e for each t e [0,T]}. 

It may be shown that the space (T>, d) is separable but not complete, however, it is possible to introduce 
a metric do on T>, equivalent to d, such that the space (T>, do) is both separable and complete. Obviously, 
we have that 

x n — > x in (V,do) <^ there exists A„ € A, sup(||x„(A„(-)) — x(-)\\e, \X n (t) — t\) — > 0. (5.5) 

The conditions imposed on the W in Theorem 14 . 1 1 invoke that E should be chosen as a separable Banach 
space continuously embedded into W '°°(r2; R 3 ). Thus we may take E = W k ' p n Wq (CI, R 3 ), where 
p € (1, oo) and k > 1 is such that 

fc>- + l. (5.6) 
P 

Thus we report the following result: 

Theorem 5.3 Let fl C R 3 be a bounded Lipschitz domain. Assume that p = p(p) satisfies hypothesis 
(|2.ip and consider the Navier-Stokes system (|1.8|) - (|1.13p where w is an E-valued Levy process, with 

E = w k ' p nw Q 1 ' 1 {n,R 3 ), fe>- + i, k> i, P > i. 

Lf the initial datum (po, (pu)a) is anXD-valued random variable then there exists a random variable (p,u) 
satisfying the Navier-Stokes system. 

The proof follows directly from Theorem 14.11 if we choose W = T> = 2?([0, T], E). It is easy to check 
that the conditions (|4.ip . (|4.2p are satisfied. 
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